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ABSTRACT 


This thesis analyzes the dynamic stability of submersible vehicles in motions 
in six degrees of freedom. A continuation algorithm is used in order to obtain the 
steady state solutions in terms of dive plane angle, rudder angle, and longitudinal 
separation of centers of gravity/ouoyancy. The equations of motion are then 
linearized in the vicinity of the above stated nominal point. The eigenvalues of the 
linearized system indicate the degree of stability of the nominal motion. The 
results demonstrate the stabilizing or destabilizing effects of general three 
dimensional motions as opposed to the traditional use of straight line level flight 
paths. Recommendations for robust design of control laws for commanded paths 


in combined horizontal/vertical planes are provided. 
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Le INTRODUCTION 


A. GENERAL 

As the versatility of underwater unmanned submersibles are 
becoming imagined and realized, the question of dynamic 
stability must be addressed and justified. Salvage, survey, 
and rescue operations as well as defense related operations of 
the future may become increasingly dependent on the ability of 
reliable unmanned underwater submersibles to perform 
prescribed missions. Theoretically, an unmanned submersible 
vehicle would not be as physically restricted in its motion as 
a manned vehicle. Because of its inhabitants and cargo, a 
manned submersible is limited mostly to straight line motion. 
This thesis seeks to explore the dynamic stability of 
submersible vehicles in motions in six degrees of freedom. 

The tight maneuvers that can be performed by an unmanned 
submersible vehicle make analysis of dynamic stability very 
difficult. The nonlinearity of the dynamic properties is very 
complicated and sometimes very hard to visualize. The most 
effective analysis of the dynamics and stability of a 
waterborne vehicle comes from observing and/or modelling the 
hydrodynamic forces acting upon it through the use of the 
translational and rotational equations of motion in six 


degrees of freedom. These equations are linearized in the 


vicinity of a nominal flight path. The eigenvalues of the 
linearized set of equations are an indicator of the relative 
degree of dynamic stability for the vehicle’s motion. 
Given a mathematical model with an established set of 
physical characteristics, analysis of the dynamic stability 
can be studied varying the dive plane angle, the rudder angle, 
and the location of the center of gravity. Analysis of the 
linearized equations of motion and determination of the 
eigenvalues by individual calculation would prove tedious. 
Therefore, computer algorithms developed by Seydel [Ref.1] and 
Aydin [Ref.2] are used in this thesis to accomplish these 
tasks. An explanation of uncoupled stability is provided in 
Chapter II. Formulation of the linearized equations of motion 
in six degrees of freedom is presented in Chapter III. 
Chapters II and III are based and partially paraphrased from 
Papoulias [Ref.3]. The experimental results in Chapter IV are 
based on the understanding of Chapters II and III. Chapter V 
summarizes the results and provides recommendations for future 


submersible modelling research. 


B. PARAMETERS 

Several variables were used in the algorithms associated 
with this research. The shorthand notation is from Smith, 
Crane and Summey [Ref.4] and is as follows: 

Variables 


Kr Va S Distances along the body fixed axes. 


u, V, W Translational velocity components of model 
relative to fluid along body axes. 


by, G #F Rotational velocity components of model along 
body axes. 


Me Vis Zi Hydrodynamic force components along body axes. 

K, M, N Hydrodynamic moment components along body 
axes. 

yw, 8, 0 Yaw, pitch, and roll angles; positive values 


following the right hand rule. 
Ray Mey Se Center of gravity coordinates along body axes. 


Man View (25 Center of buoyancy coordinates along body 
axes. 


vr Iz, Moments of inertia about body axes. 

The standard body-fixed, right-hand orthogonal axis system is 
employed throughout all data simulations. All data run 
simulations described herein were initialized with a constant 
propeller speed of 500 rpm, which translates to a value of 6.0 
ft/sec for "u". The data is analyzed for rudder angles 
between 0 and 20 degrees and stern diving plane angles between 
-20 and 20 degrees. Separate data runs were conducted within 
the above stated parameters for simulated model configurations 
with the center of gravity shifted longitudinally forward and 


aft between -1.5 and 1.5 percent of the model length. 


II. UNCOUPLED STABILITY 


A. MANEUVERING IN THE HORIZONTAL PLANE 
1. Basic Maneuvering 

It is extremely important to note that the ability to 
control a vessel involves two major characteristics that are 
contrary to one another. These two characteristics are the 
ability of controls fixed course stability and the ability of 
the vessel to turn. Obviously a vessel that is extremely 
course stable will not turn easily. Newton’s equations of 
motion, used with reference frames both fixed with respect to 
the earth (inertial) and fixed relative to the vessel will 
illustrate the dynamics of ship motion. This motion will be 
first described in the horizontal plane and later described in 
the vertical plane. As a body moves within the inertial 


frame, its motions can be described as follows: 


m xX, = X, (1) 
mY, = Y, (2) 
I,¥=N (3) 


through the use of the variables introduced in Chapter I. The 


subscript o pertains to motion in the inertial reference 


frame. In the vessel’s relative frame of reference, these 


equations become: 


xX =X, cosy + ¥, siny , (4) 
Y= Y, cos - X,siny , (5) 
X, = ucosy -vsiny, (6) 
Yo = usinp + vcosy, (7) 
and 
ber. (8) 


Therefore, equations (6) and (7) can be integrated and 
substituted into equations (1) through (3). They in turn can 
be substituted into equations (4) and (5) with the following 


results: 


mu-mvr =X, surge equation, (9) 
mv +mur=¥Y, sway equation, (10) 
I,f=N, yaw equation. (2) 


These equations are true only if the reference plane center 
and the vessels center of gravity coincide. If the two are 


indeed located at different locations, the following equations 


result: 
mu ~ mvr - mx,r2 =X, (12) 
mv + mur + mx,r =Y, (13) 
I,t + mx,(v + ur) =N, (14) 
provided that y, = 0. An important point to note is that 


centrifugal forces do not appear in equations (1) through (3). 
These forces are not evident until the relative reference 
plane is used in equations (9) through (12). Centrifugal 
forces are not experienced on the inertial reference frame 
Since it is stationary but are indeed evident when the 
relative reference plane of the vessel is studied. The 
hydrodynamic forces X and Y and the moment N are composed of 


four different types of forces: 


1. F forces; Fluid forces acting on the hull by the 


surrounding fluid (water). 


2. R forces; Control surface forces due to rudders, bow 


planes, dive planes, thruster, etc. 


3. E forces; Environmental forces due to wind, current 
and waves. 
4. T forces; Propulsion forces due to propellers, 


thrusters etc. 


Due to the construction characteristics of most waterborne 


vessels today, the following can be stated: 


X=Xp,+ Xpt+X,p+T, (15) 
Y=ryYet+ Yet Vp « (16) 
N=N,+ Na+ Ne . (17) 


Since the hydrodynamic forces X,, Y,;, and N,; depend on vessel 
motion and must therefore be functions of vessel velocity and 


acceleration through the water we can write: 


Xp t Xpluevit, Vv, 0,4) , (18) 
Y, = Y,(u,v,u,V,7,£) , (19) 
Nz = Np (u,v,uU,v,r,2) . (20) 


The interrelationship of equations (18) through (20) is a 
complicated one. Usual maneuvering studies simplify this 
interrelationship. The object is to study ship response about 
a nominal equilibrium point (designated by the subscript 1). 
Through the use of the Taylor series expansion around the 
nominal point, and keeping only the first order terms, the Y 


relationship becomes: 





Be oi ; Y, oY, 
Y, = Ye(U,,Vye 0s VyeFy02,) + (u- u,) at + (v - Mi) ae 
oy, 
# eae (=> YS , (21) 
VOY; 


where all of the partial derivatives are evaluated at a 
nominal set of conditions. X, and N, can be similarly 
manipulated and evaluated. 

The nominal set of conditions is the set of variables 


associated with straight line motion: 


u, =U, v, =u, =v, = 27, = 24, =0 (22) 


For the sake of understanding, the vessel is considered to be 


port/starboard symmetric. Therefore 6F/déu = 6F/d6t0 = 0 since 


a change in the forward velocity or acceleration will produce 
no transverse force in vessels that are also symmetric in xz 
plane. In addition, if the vessel is indeed in equilibrium, 
and in straight line motion, there can be no Y,-forces. 
Therefore the first term in equation (21) is also zero. Once 
these additional considerations are made, equation (21) 


reduces to: 


by, Oe. Oe. 8%, 


= —_ = TE 23 
oO ae ae’ aE ee ea 


In like fashion, the surge force and the yaw moment can be 


simplified: 
OX, OX, OX, OX, OX, Ox, 
= a - aa peal pala ak 24 
xX, att (u - U)+ ae Fa oT tr, (24) 


Where the cross coupling derivatives 6Y,/6r, 6SY,/dZ, ON,/dv, 
and 6N,/6v¥ usually have small nonzero values because most 
vessels are not symmetrical about bhi 4m plane even if that 
plane is at the midlength of the vessel (bow and stern 


contours are usually drastically different). However the 


cross coupling derivatives 6X,/6r, 6X,/6f, 6X,/év, and 6X,/6v 


are zero because of the symmetry about the xz plane and 


equation (22). Equation (24) now reduces to: 
Ox 
Xp = —u + OR i -U) . (26) 
u u 
Using standard notation, 6Y,/év = Y,, O6N,/6r = N, ... etc. 


These are known as the hydrodynamic derivatives. The physical 
meaning of these mathematical derivatives is very important to 
the development of understanding the connection between 
maneuvering and hydrodynamic forces. For example, Y, is the 
force in the sway direction due to a unit change in the yaw 
angular velocity. 

Using this conventional notation, and substituting into 
equations (12), (13), and (14), the linear equations of motion 
in the horizontal plane in the absence of environmental 
disturbances (ideal handling conditions) and with the control 


surfaces set at zero, become: 


(m- X,)u=X,(u-0) , (27) 


(m- ¥,)v - (¥, -mx,) $= Yv+(y, -mu)r , (28) 
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(I, - Nz) ft - (Ny - mx,) Vv = Nv+ (N 


E 


- mx,U)r , (29) 


where we can see that within linearity the surge equation 
decouples from sway and yaw. 

Note that all of the terms of equations (28) and (29) must 
include the effect caused by holding the vessel’s rudder at 
zero in order for all the previously imposed simplifying 
conditions to hold true. If the possibility of non-controls 
fixed motion were to be considered, equations (28) and (29) 
must include terms on their respective right sides expressing 
the forces and moments caused by control surface deflection 
(i.e. rudder deflection) as a function of time. Assuming that 
the control surface most generally used in the horizontal 
plane is indeed a rudder and that the rudder force and moment 
on the vessel are functions of the rudder angle 6 only, and 
not §, the following is true: 


OXn 5 


Xp(5) = X,(8=0) + 38 


= 69 a (30) 


o) + Reeves , (31) 


¥_(3) = 


YR (6 


an, 
0) + S86 = ms . (32) 


_ Ng(8) = Ng lS 


14. 


Positive rudder deflection corresponds to a port turn when 
associated with a stern rudder. Ys and Ns are the rudder 
hydrodynamic derivatives. If the rudder forces and moments 
are now included in the linearized equations for sway and yaw 
(equations (28) and (29), respectively), the equations of 


motion become: 


(m- Y,)v - (¥, - mx.) f= Yyv+ (¥, - mU)r+ ¥,5 , (33) 


(I, - Ny) rf - (N, - mx,)v = N,v + (N, - mx,U)r + Ngo. (34) 
2. Stability of Straight Line Motion 

The concept of an object remaining on a set path is 
related to the concept of a vessel’s course stability or the 
stability of its straight line motion. There are various 
types of motion stability associated with marine vehicles and 
they are separated by the attributes of their initial 
equilibrium state that are retained in the final path of their 
center of gravity. These types of motion stability are 
depicted in Figure 1. For each type of stability depicted, 
the vessel is initially travelling at constant speed along a 
straight path. Case I, straight line or dynamic stability, 


retains the attribute of straight line motion but not the 
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CASE | 
STRAIGHT-LINE STABILITY 


CASE IIA 
OIRECTIONAL STABILITY 
LOW DAMPING 


CASE 118 
OIRECTIONAL STABILITY 
HIGH DAMPING 





CASE I! 
POSITIONAL STABILITY 








— — — — — INTENDED PATH 





ACTUAL PATH 


pes POINT OF DISTURBANCE 


Figure 1: Various Kinds of Motion Stability 


— 
WwW 


attribute of direction from the initial state of equilibrium 
after release from a disturbance. In Case II, directional 
stability, the final path after release from a disturbance 
retains not only the straight line attribute of the initial 
path, but also its direction. Case III is similar to Case II 
except that the vessel does not oscillate after the 
disturbance, but passes smoothly to the same path as in Case 
Tis Finally, in Case IV, positional motion stability, the 
vessel returns to the original path and retains direction as 
well as position. 

The types of stability are classified in an ascending 
order. A directionally stable vessel must also posses 
straight line stability. One that is positionally stable must 
posses both directional and straight line stability. Straight 
line stability or instability is indicated by the solution to 
a second order differential equation, directional stability or 
instability is indicated by the solution to a third order 
differential equation, and positional stability or instability 
is indicated by the solution to a fourth order differential 
equation. 

Each type of controls fixed stability has an associated 
numerical index which by its sign determines the relative 
stability or instability of the ship. The magnitude of the 
index determines the degree of stability or instability. “In 
order to illustrate the determination of this index, and to be 


able to evaluate the degree of dynamic stability, a return to 
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the differential equations of motion (28) and (29) is 
necessary. This homogeneous set of differential equations 


will have a solution of the form: 


v(t) = yer + yer , (35) 
r(t) = re" + re , (36) 


where the constants vi, V2, Y;, XY, are determined from the 
initial conditions. Differentiation and subsequent 
substitution of equations (36) and (37) into equations (28) 


and (29) yields the following characteristic equation: 


Ao +Bo+C=0 , (37) 


where 
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A= (I, - N,) (m- ¥,) - (Y¥, - mx,) (Ny - mxX4) , (38) 
B= -(I,- N,)¥,- (m- Y,) (N, - mx,U) 
- (Y¥, - mU) (Ny - mx,) - (Y¥, - mx.) NM, , (39) 
C= (N, - mx,U)Y,- (Y¥,- mUu)N, , (40) 
with roots 
- 1,-2 By2 _ gil (41) 
01,9, a ar (3) 4(5) Jog 


If on course stability is to be displayed, then o, and o, will 
have negative real parts. Referring back to equations (36) 
and (37), both v(t) and r(t) will go to zero for increasing 
time. This means that the ship will eventually assume a 
straight line direction at a different, in general heading. 
This is straight line stability as defined in Figure 1. 


Symbolically, the stability requirement 


R{o,,0,} < 0, 


translates to 
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with all the coefficients of the quadratic equation (38) are 
positive. To evaluate these conditions, we must check the 
signs of the hydrodynamic derivatives that make up A,B and C. 
A physical interpretation of the hydrodynamic variables as 
forces and moments acting on the hull due to the various hull 
motions with respect to the water is necessary. 


After some analysis, the following conclusions are drawn: 


always negative and large, 
always negative and large, 
always negative and large, 
always negative and large, 
and of uncertain sign, 
and of uncertain sign, 
and of uncertain sign, 


and of uncertain sign, 





small and of uncertain sign, 
always positive and large, 


always positive and large. 


With these conditions described thusly, A and B are therefore 


always positive and stability must require: 





C= (N, - mx,U)¥, - (¥, - mUu)N,>0 . (42) 
This is the Stability Criterion. The more positive C is, the 
more stable the ship will become and the more increasingly 
difficult the ship will be to maneuver. The more negative C 
is, the more unstable the ship will become and the continuous 
- use of rudder will be required to maintain course. Although 
X, 1s small and of uncertain sign, for stability it is desired 
to have x, > 0 (center of gravity forward of amidships). 


Standard ship design experience agrees with this observation. 


B. MANEUVERING IN THE VERTICAL PLANE 

1. Linear Equations in the Vertical Plane 

Maneuvers in the vertical plane can be describe in much 
the same manner as they can in the horizontal plane. The 
logical progression in the determination of stability will be 
traced out here for the vertical plane as it was in the 
previous section for the horizontal plane. 

A submerged vessel moving in the vertical plane has a 
nominal path defined as the level path at the commanded depth. 
It has a pitch angle, 8, which is defined as the angle between 
the vehicle’s longitudinal axis and the x axis and an angle of 
attack, @, which is defined as the angle between the x axis 
and the total velocity vector. Figure 2 illustrates these 


points. If Newton’s law of motion is once again used, the 
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Nw =< 





Figure 2: Vehicle Geometry in the Vertical Plane 


following equations are true using the variabies introduced in 


Chapter I: 


mt, =X , (43) 
m3, = 2 , (44) 
roe, (45) 


As in the horizontal plane, the vertical transformation of the 


coordinate system yields: 
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X, = +Xcos@ + Zsin6é , 
Z, = -XsinO + Zcos6 , (46) 
X, = +ucos6 + wsiné , 


Z = -usin6 + wcosé , 


where the equations of motion are: 


mu+mwq=xX , surge equation , 


mw ~ muq=Z , heave equation , (47) 


Ija=M , pitch equation , 


and 


@=q. (48) 


When expressed with respect to the vessel’s reference frame, 


these equations become: 


mu + mwq ~ mx,g? + mx, =X , , (49) 
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mw - muq - MxX,q - mx,q? =Z , (50) 


Iq - mx,(w - ug) + mx,(U- wq) =M , (51) 


where (x,,Z,) are the coordinates of the ship’s center of 
gravity with respect to the ship fixed reference frame, and we 
assume y, = 0. 

Equations (49) through (51) demonstrate a coupling between 
surge motion and vertical plane motion that was not present in 
the discussion on the horizontal plane. If the vertical plane 
equations of motion are linearized around a level flight path 
at constant speed, U, we can see from equations (49) and (51) 
that the two terms mxq and mx will remain after 
linearization. This means that heave and pitch are coupled to 
surge even for small motions. This dynamic coupling arises 
due to the nonzero z, term. When a vessel pitches, a nonzero 
surge acceleration will contribute to the inertial moment. 
Obviously there was no such linearized dynamic coupling in the 
horizontal plane. Although vertical plane motion must be 
studied together with surge, the coupling discussed here is 
relatively small due to small values of z, and in this case 
can be neglected. Therefore, the uncoupled linearized 


equations of motion in the vertical plane are: 


mw- mUq-mx,q=Z , (52) 
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Iq - mx,(w- Ug) =M . (53) 


Much like in the previous section, Z and M can be expressed in 
terms of slow motion derivatives. Using vertical plane 


notation they are: 


Z=Zw+Zoq+Zw+ Zq+Z,8 , (54) 


M=Mw+Mq+M w+ Mqd- (z,- Z,;)Wsind+Mé , (55) 


where Z, stands for the change in the hydrodynamic force in 2 
with a unit change in w (note: Z, is synonymous to Y, in the 
horizontal plane) and the same for the other variables. The 
additional term in equation (55) is the hydrodynamic restoring 
moment, and all other variables are in accordance with the 
listing given in Chapter I. Provided the vessel is properly 
balanced (weight = buoyancy), there is no such hydrostatic 
force in the heave equation (54). The 6 in these equations 
refers to the stern plane angular deflection as shown in 
Figure 2. 

Equations (52) through (55) lead to the linear uncoupled 


equations of motion in the vertical plane: 


(m- Z,)w- (mx, + Zy)q = Zw + (Z, + mU)q + 2,8 ; (56) 
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(I, - My) @ - (mx, + My) w= Mw + (M, - mx,U)q - (2, - 23) W6 + 4d 
(57) 


These two equations can be used with equation (48) in order to 
predict the response of the vessel in the vicinity of a 
nominal level flight path. Deviations in depth can be 


computed from: 


z = -Usin® + wcosé , (58) 


or 


z= -U06 (59) 


in its linearized form. 
2. Stability of Motion 
Stability of motion in the vertical plane can be 
studied by using equations (56) and (57) after 6 is set to 
zero. As in the case of the horizontal plane, small motion 
are being studied around the nominal flight path with the 


controls fixed. Equations (56) and (57) become: 
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(m - Z,)w - (mx, + Zy)q - Zw - (m+ Z,)q = 0 , (60) 


(Ly - Mg) @ -(mx, + M,) W -M,w -(M, - mxg)q- MO=0 , (61) 


where 


My = -(Z - Zg)W=-(2,- Z,)mg , (62) 


and U = 1 since dimensionless units are commonly used. The Mg 
term is the hydrostatic "spring constant" due to the nonzero 
metacentric height. The corresponding term in the horizontal 
plane is M, and is equal to zero. The negative sign in 
equation (62) makes M, < 0, which is consistent with 
established definitions for positive forces and moments in the 
right handed orthogonal coordinate system. Manipulating 
equations (60), (61), and (48) and transforming to the s 


domain gives the following: 


((m - Z,)s - Z,]w +[-(Z, + mx,)s? - (Z,+m)s]6 , (63) 


[-(M, + mx,)s - M,]w+[(I, - My) s? -(M, - mx,)s -M)]@=0 , 
(64) 


where the characteristic equation is: 
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As? + Bs*+cCs+D=0 , (65) 


where, 
A= (m- Z,) (I, - My) - (My, + mx.) (Zy + mxX_) , (66) 
B= -(m- Z,) (My - mxg) - (Zy + mx) M, - (Z, + m) (M, + mx,) 
(I, -M)Zy s (67) 
C= (My - mx,)Z, - (24g + mM, - (m- Z,)Mm , (68) 
D=Z,My . (69) 


When comparing equations (66) through (69) with equations (38) 
through (40), coefficients A through C are similar in the 
horizontal and vertical planes with the extra metacentric 
moment term appearing in equations (68) and (69). If there 
were no hydrostatic terms, this cubic characteristic equation 
(65) would reduce to a quadratic similar to the equation 
derived in the horizontal plane. This means that there is a 


zero root in this case. A physical analysis is necessary to 
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understand the full meaning of this mathematical fact. In the 


vertical plane: 


(As? + Bs? +Cs+D)O=0 . (70) 


This means that, assuming certain conditions of stability are 
satisfied, it is possible for 0(t) ~0 as to. Therefore, 
directional stability is possible. Where in the horizontal 


plane, 


(As? + Bs + C)sp = (As? +Bs+c)f=0 . (71) 


This means that the best that can be expected in the 
horizontal plane is W(t) = r(t) ~0 as t 4, which means that 
yw will be increasing linearly with time. Therefore, 
directional stability is not possible and the best possible 
situation is straight line stability. 

In the case of the cubic equation the Routh-Huritz 
conditions for stability dictate that all coefficients A, B, 


C, and D be positive and that: 
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BC-AD>O . (72) 


After close analysis similar to that conducted in the 
horizontal plane, it can be proven that A, B and D are 
positive. Therefore in order for the vessel to stable in the 


vertical plane: 


C= (My - mx,)Z, -(Z, + m)M, -(m - Z,) My > 0 


Since (m - Z,) > 0 and Mg < 0, this condition is easier to 
satisfy in the vertical plane than in the horizontal plane. 
This proves that a positive metacentric height contributes to 
dynamic stability as well as static stability. Assuming that 
equation (73) is satisfied, it can be shown that equation (72) 
is always satisfied. This is best shown by using a method 


detailed in [Ref.3]. If the coefficients are written as: 


where the subscript 0 refers to the hydrodynamic coefficients 


without’ the hydrostatic terms, and subscript 1 refers to the 





metacentric terms. The parameter € is the metacentric height 


(Zo - Zp). Therefore: 


BC - AD = BC, + €(ByC, - AyD,) . (74) 


Based on the given conditions, By, C,, Aj), and D, are all 
positive. Therefore the difference (B,C, - AD,) is a small 
number of uncertain sign. Since it is multiplied by €, a 
small positive number, the result is a small number of 
uncertain sign. This small number, regardless of sign will not 
change the sign in equation (74), assuming that C is 
positive. This means that Cy > O will be a sufficient 
condition to satisfy vertical plane stability. The vertical 


plane stability reduces to: 


(My - mX_)Z, —(Zg + m)M,> 0 . (75) 


Where this is a sufficient condition for vertical stability, 
it is not a necessary condition. It is possible that the 
hydrodynamic term will stabilize the vessel and equation (72) 
will be satisfied (and all coefficients will be positive), 
even though equation (75) might not hold true. This condition 


is dependent on the magnitude of the metacentric height and 
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the speed of the vessel. As a vessel experiences an increase 


in forward speed, it will behave less stably in the vertical 


plane. 

Recalling the Stability Criterion for the horizontal 
plane from equation (42), equation (75) can serve as the 
Stability Criterion for the vertical plane. These two 


expressions are commonly normalized into Stability Indices: 
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G, 1 Z,, (Mg MXq) ’ ( ) 
N,(Y, - m) 

G, = ees 77 

" Y,(N, - mx) ay 


These two indices are the vertical and horizontal plane 
stability indices and they are used to asses maneuverability 
and relative stability. They are always less than one, and 
positive values indicate stability and negative values 


indicate instability. 
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IIIT. MANEUVERING IN SIX DEGREES OF FREEDOM 


A. RIGID BODY EQUATIONS OF MOTION AND KINEMATICS 
Development of the general maneuvering equations in six 
degrees of freedom requires the employment of two coordinate 


systems depicted in Figure 3: 


OC, YZ). = reference inertial frame, 
AUK ¥ y's) % moving frame fixed on the vessel’s 


geometric center. 


(0) 
xX 
| . = movin 4 
Pp VY |t 
y 2 
4 cE 


Figure 3: Fixed and Moving Coordinate Systems 
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As discussed in the previous chapter, Newton’s law requires an 
inertial frame, while the motion of the vessel is best 
described in its own frame. G denotes the center of gravity 
of the moving body, which is lcoated at (x,,y.,,Z,) with respect 
to the A(x,y,z) frame. Unit vectors in this moving frame are 
denoted by 4,f,8, in the x,y,z directions, respectively. 

The translational equations of motion can be found through 
the use of the principle of linear momentum in the inertial 
frame. The rate of change of the linear momentum equals the 
sum of the external forces. Vector analysis is then used to 
express these quantities in the A system. Chapter 4.1 of 
Reference 3 explains this conversion in detail. Once the 
conversion to the A system is accomplished, Newton’s law can 


be broken into components as: 


m(u + qw- rv - X,(q* + ©?) + yg(pq@- ft) + Z26(pr+q@)] =x , 
(78) 

m(v + ru - pw- y,(r? + p*?) + zg(ar - Bb) +x,(apot+ ry] =Y, 
(79) 

mlw + pv - qu- Z,.(p? + q*) + x,(rp - @) +yglrq+pb)] =Z , 
(80) 


where the total exciting force ,F, is denoted by, 


ent 


Fa xPe vyfu ge . (81) 


Therefore the external moment,M can be written as: 


M= KE+MP+NK , (82) 


and broken into components as: 


I, p+ (Ll, > £)@r + £,(pr = @) = £,(¢q* = 27) = (pq + £) 
+ mlyg(w - uq+ vp) - z.(v¥ + ur - wpo)] =K , ~~ (83) 
22+ (he, - Tip - flee +l + hie 2h + 2, te - 24d 
-m[x,(w - uq + vp) - z,(u- vr +wq)] =M , (84) 


beh + Sy Dee - 2 is" + gy: = Ling + a) + hee - 


+m[x,(v + ur - wp) - yg(t- vr +wq)] =N . (85) 


Equations (78) through (80) are the translational equations of 
motion, and together with the rotational equations (83) 
through (85), they make up the six degrees of freedom of 


motion of a rigid body expressed in a coordinate system moving 
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along with the body. Once again the variables involved are 
the same variables introduced in Chapter I. 

Although the equations of motion are written in a frame 
associated with the vehicle, the vessel’s orientation and 
position in the inertial frame is still very important. The 
angular (orientational relationship between the two frames is 


determined by the three Euler angles: 


yaw or heading,wy : rotation about x axis 
pitch or elevation,@ rotation about y axis 
roll or bank, : rotation about z axis 


Any vector described in the two reference frames are related 
by some transformation matrix, T which is a function of these 


three Euler angles. 


Therefore, 


T=T(o,0,4) . (86) 


The general transformation matrix can be broken up into three 
individual rotations. The order of application is very 


critical in order to maintain the integrity of the 
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transformation. Physically, infinite rotations do not 
commute. The prescribed manner is to first yaw, then pitch, 
and finally to roll. Standard matrix notation therefore 


yields: 


T(o>,6,y) = T(h) T(B) Tip). (87) 


The full transformation matrix is: 


cos6cosy cosOsiny -sin@é 
T(>,8,) =|-sinycosd+sindsinOcosy cosdcosy+sindsinésiny sindcos6 |(88) 
sindsiny+cosdcosysiné -sindcosy+cosdsinOsiny cosdcos6 


Small rotations do commute, and for small angles, 


1 yw -6 
T(>,6,y) =|-w 1 oo] . (89) 
6 -) 1 


and in this case T(6,89,y) = T(8,6,W). Equations (87) and (88) 


can now be used to transform any vector in one system to its 
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equivalent in the other system through the use of matrix 
multiplication. Here are some important examples of the 


relationship between the two reference frames: 


X, = ucos@Ocosy + v(-sinycosd + sindsinfcosy) 
+w(sin6siny + cos@cosysiné) , (90) 


Yo = ucosOsiny + v(cosdcosy + sindsinOsiny) 


+w(-singdcosy + cosdsin@siny) , (91) 

Z, = -u sin® + v sindcos@ + wcosdgcosé , (92) 
p=-sinO+o , (93) 

gq= sindcos@ +O cosh , (94) 
r=wcos@cosd - 6 sind , (95) 
@=p+qsindtan6 + rcosdtané , (96) 

6 =qcoshd - rsind , (97) 

w = ging , cos (98) 


Where the subscript o represents the inertial reference 


frame. For small angles, equations (96) through (98) become: 
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which means that the rates of change of the Euler angles are 
the same as the angular velocities only for infinitesimal 
rotations. The six kinematic equations (90) through (92) and 
(96) through (98) along with equations (78) through (80) and 
(83) through (85) provide a complete description of the body’s 


motion in six degrees of freedom. 


B. EXTERNAL FORCES AND MOMENTS 

Equations (78) through (80) and (83) through (85) are the 
six degrees of freedom equations of motion for a rigid body 
where the left hand terms represent inertial forces and 
moments while the right hand terms represent the external 
forces and moments. These external forces and moments include 
hydrodynamic contributions as well as weight and buoyancy 
effects and forces from control surface deflections and 
propeller thrust. Hydrodynamic restoring forces and moments 
are due to the vehicle weight and buoyancy. Forces and 
moments due to control surface deflections are experienced as 
added drag in surge, while in sway, heave, pitch, and yaw they 
are directly proportional to the deflection. Just as the 


position and velocity vectors in the previous section can be 
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expressed in terms of either the inertial reference frame or 
the vessel’s relative reference frame, so can these forces and 
moments. The tie between the two frames of reference in, once 
again the transformation matrix. Particularly in the case of 
the hydrodynamic force and moment contributions, the 
derivation of these binding expressions can be very 
complicated and lengthy. If a more in depth discussion of the 
subject is desired, Reference 3 provides an excellent 


overview. 


Cc. LINEARIZATION 

Linearization, as taught as a mathematical tool is useful 
for understanding the behavior of nonlinear systems. The 
principle is used here to analyze particular vehicle nonlinear 
motion about a nominal flight path. Generally, linearization 
is performed in the vicinity of a nominal point. This point 
is that point where the system is expected to spend most of 
its life. The nominal point is defined by the system 
condition that all time derivatives equal zero. Any existing 
system controls can be fixed or expressed as a function of the 
nominal point. Once a system solution has been computed using 
the nominal point and the necessary qualifing condition, the 
system can be linearized. This is done by determining the 
Jacobian matrices of the system. This new linearized system 
can be used to analyze system open loop stability through the 


examination of the eigenvalues of the Jacobian. 
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Iv. DISCUSSION AND RESULTS 


A. GENERAL 

The major motivation of this thesis is to further explore 
the line of reasoning that underwater vehicles are no longer 
bound by the straight line motion that is caused by personnel 
or cargo staged onboard the vessel. This thesis seeks to 
prove/verify that a particular underwater vehicle is 
stable/unstable when enough physical information is known. 
While gathering information to form the stability analysis, 
valuable information can be obtained concerning the motion of 
the vehicle while it performs three dimensional maneuvers. 
The objective of this study is to provide a parametric 
stability analysis of the Swimmer Delivery Vehicle in six 
degree of freedom motion. All physical .properties of the 
vehicle are known as well as the linear, nonlinear and 
coupling terms. The main object is to asses control design 
while keeping the controls fixed and varying the distance 
between the longitudinal center of gravity (LCG) and the 


longitudinal of buoyancy (LCB). 


B. PROCEDURE 


After a prerequisite understanding of uncoupled motion in 
the horizontal and vertical plane, the general procedure 


followed in this parametric study continued with the full six 
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degree of freedom equations of motion. The necessary 
variables in order to analyze the system dynamics are the 
translational velocities surge, sway, and heave (u, v, w); 
rotational velocities roll, pitch, and yaw (p, qd, r); the 
Euler angles of roll and pitch (o, 8); and the rate of change 
of heading - the first derivative of the yaw angle, Y. In 
order to provide a comprehensive analysis of vehicle behavior 
and capability, a wide variety of realistic control parameters 
and loading conditions were chosen and evaluated for their 
contribution to the steady state motions of the above stated 
variables. The final step in the procedure is to linearize in 
the vicinity of the steady state and to compute the system 
eigenvalues. These eigenvalues serve as indicators of the 


relative stability of the vehicle under the given conditions. 


C. STEADY STATE ANALYSIS 

The steady state analysis started with the full, 
nonlinear, coupled six degree of freedom equations of motion. 
All time derivatives were set to zero since steady state is 
defined as that point where there is no change in the state of 
a variable with respect to time. The resulting algebraic 
system of coupled, nonlinear equations was solved for the 
control parameters and loading conditions used. In order to 
perform an adequate analysis, a continuation algorithm 


developed by Seydel [Ref.1], with subroutines by Aydin [Ref.2] 
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were used. The following range of control parameters were 


examined: 


- stern diving plane angle varied between +20 
degrees 

- rudder angle varied between 0 and 20 degrees 

- LCG/LCB separation between +1.5 percent of vehicle 


length 


The graphical results are organized in sets according to 
LCG/LCB separation. The LCG/LCB separation was varied between 
+1.5 percent of vehicle length in increments of 0.5 percent. 
Each set of the following continuation results depicts the 
entire stern diving plane angle range as the x axis. Each set 
has a plot where the y axis is one of the aforementioned 
variables u, v, w, p, q, xr, 6, 8, and YW. The translational 
velocities are plotted in ft/sec, and the angular velocities 
in rad/sec. All simulations were run at a constant propeller 
speed of 500 rpm, which translates to a value of 6.0 ft/sec 
for the forward speed, u. The final plot in each set is the 
relative degree of stability which will be discussed in the 
next section. These plots must be analyzed as a set in order 
to provide a comprehensive understanding of the vehicle six 
degree of freedom motion. Through the comparison of these 
plots, the motion associated with any maneuver within the 


range of parameters can be predicted. 
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De STABILITY ANALYSIS 

Steady state analysis of the maneuver alone is not 
sufficient. It must be determined if the vehicle is to remain 
stable (controllable) while the maneuver is negotiated. 
Therefore the final plot of the set depicts the degree of 
stability plotted against the range of the stern diving plane 
angle. The units of the degree of stability are in terms of 
inverse time, sec. This index is obtained by evaluating the 
9 x 9 Jacobian matrix created by the determination of the 
steady states discussed in the previous section. The largest 
real part of all resulting eigenvalues in a continuation is 
the degree of stability. This, in physical terms, is the most 
dominant time constant for the system. The index measures the 
slowest exponential convergence to the steady state motion 
when negative or the fastest exponential divergence from the 
steady state motion when positive. The more negative the 
index, the more stable the vehicle maneuver. Conversely, the 


more positive the index, the more unstable the maneuver. 


E. GRAPHICAL RESULTS 

The following figures along with an understanding of the 
previous sections of this chapter illustrate to the reader the 
dynamics of vehicle motion and the stability of that motion 
within the parameters inputted into the continuation 


algorithm. 
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Figure 5: v vs Stern Plane Angle; x, = 0% 
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Figure 6: w vs Stern Plane Angle; x,,= 0% 
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Figure 7: p vs Stern Plane Angle; x,,= 0% 
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Figure 9: r vs Stern Plane Angle; x,g= 0% 
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Figure 11: 8 vs Stern Plane Angle; x, = 0% 
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Figure 13: Stability Index vs Stern Plane Angle; xX,,= 0% 
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Figure 15: v vs Stern Plane Angle; 
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Figure 17: 


p vs Stern Plane Angle; x,,=.5% 
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Figure 21: @ vs Stern Plane Angle; 
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Figure 22: W vs Stern Plane Angle; x,,=.5% 
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Figure 23: Stability Index vs Stern Plane Angle; 
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Figure 24: u vs Stern Plane Angle; x,,=1.0% 











Figure 25: v vs Stern Plane Angle; x,,=1.0% 
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Figure 26: w vs Stern Plane Angle; x,,=1.0% 
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Figure 27: p vs Stern Plane Angle; x,,=1.0% 
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Figure 28: q vs Stern Plane Angle; X,,—=1.0% 
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Figure 29: r vs Stern Plane Angle; x,,=1.U% 
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@ vs Stern Plane Angle; x,,=1.0% 
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Figure 32: W vs Stern Plane Angle; x,,=1.0% 





Figure 33: Stability Index vs Stern Plane Angle; X,g,=1.0% 
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Figure 34: u vs Stern Plane Angle; x,,=1.5% 





Figure 35: v vs Stern Plane Angle; x,,=1.53 
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Figure 36: w vs Stern Plane Angle; xX,,=1.5% 
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Figure 37: p vs Stern Plane Angle; x, ,=1.5% 
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Figure 38: q vs Stern Plane Angle; x,,=1.5% 
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Figure 39: © vs Stern Plane Angle; x,,=1.5% 
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Figure 40: 6 vs Stern Plane Angle; x,g=1.53% 
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Figure 41: 6 vs Stern Plane Angle; x,,=1.5% 
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Figure 42: % vs Stern Plane Angle; x,g=1.5% 
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Figure 45: v vs Stern Plane Angle; x,g=-.5% 
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Figure 46: w vs Stern Plane Angle; x,,=-.5% 





Figure 47: p vs Stern Plane Angle; x,,=-.5% 
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Figure 48: q vs Stern Plane Angle; x,3=-.5% 
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Figure 49: r vs Stern Plane Angle; x, =-.5% 
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Figure 50: ® vs Stern Plane Angle; x,,=-.5% 
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Figure 51: 8 vs Stern Plane Angle; x,,=-.5% 
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Figure 52: % vs Stern Plane Angle; 
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Figure 54: u'vs Stern Plane Angle; x,,=-1.0% 
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Figure 55: v vs Stern Plane Angle; x,,=-1.0% 
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Figure 57: p vs Stern Plane Angle; x,,=-1.0% 
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Figure 60: ® vs Stern Plane Angle; x,,=-1.0% 
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Figure 61: 8 vs Stern Plane Angle; Mayge= 1.03 
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Figure 63: Stability Index vs Stern Plane Angle; x,g=-1.0% 
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Figure 64: u vs Stern Plane Angle; x,,=-1.5% 
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Figure 65: v vs Stern Plane Angle; X,g=-1.5% 


72 








ee ee ee er ee te ee eee ee 


-20 “15 -10 “5 0 5 10 15 20 











Figure 67: p vs Stern Plane Angle; xgg=-1.5% 
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Figure 69: rvs Stern Plane Angle; x gg=-1.5% 
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Figure 70: @ vs Stern Plane Angle; x,,=-1.5% 




















Figure 71: 9 vs Stern Plane Angle; x,,=-1.5% 
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Figure 72: ~ vs Stern Plane Angle; x_=-1.5% 


ci 


\ 
| 
| 


Le | 
SOS gl ee 4 
eS i 


-0.064 os 








-0.066 es 


wwe. 


' 


Saewas: een, (2) oneness Domne 1 








Figure 73: Stability Index vs Stern Plane Angle; “Xgg=- 1.5% 
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V. CONCLUSIONS AND RECOMMEDATIONS 


The results of the previous chapter show that the 
continuation algorithm is efficient in tracing the complicated 
steady state solutions associated with motions in six degrees 
of freedom. The translational and angular velocity 
components, the Euler angle components for roll and pitch and 
the and the rate of change of the vehicle’s heading are all 
products of the algorithm. The degree of stability index was 
also presented in graphical form in order to determine if the 
vehicle would remain under positive control while maneuvering. 
The more negative the index, the more stable the vehicle 
motion. In every set of data the index became more negative 
as rudder angle was increased. Therefore, for each set, the 
straight line case was the least stable. This same effect, 
though not as pronounced, is present for non zero dive plane 
angles. Recommendations for future research in related areas 
include the evaluation of different control schemes with 
respect to performance. Some examples are the verification of 
vehicle path accuracy, the ability of the vehicle to reject 
external disturbances and remain on course, and the degree of 
confidence in performance and stability if there are errors 


between the physical and mathematical modeling of the vehicle. 
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